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Abstract 

(N 

' A series of high presicion atomic experiments was carried out last decade to 

get a better estimate for the electron (anti) neutrino mass m u . The reaction 
to be observed is molecular tritium (5- decay. The m u value serves as one of 
the tune parameters to fit (3- electrons spectrum via reference theoretical one. 

o 

\q ' The unexpected message consists of that parameter has to be fixed at 

o ' 

some negative value to get best statistically reliable fit. Apart of some exotic 

I 



scenarios suggested the problem remains open. In this paper we tried to 
reanalyze ground features of the final states spectrum (FSS) and its influence 
on the (3- electrons spectrum. A new approach has been developed which 



gives some hints to make proper modifications to the tabulated FSS in the 



^ ■ course of experimental fits. 

I. INTRODUCTION 

The last atomic neutrino experiments base on the observation of molecular tritium j3- 
decay [|TJ : 

T 2 -> T 3 He + + e~ + z> e (1) 

The high energy edge of j3- electrons spectrum corresponds to about 18.6 KeV. The sudden 
change of the nuclear charge and hard recoil kick lead to the population of the hundreds 
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electronic- rovibrational states of final T 3 He + molecule. The neutrinos v e are hardly observ- 
able but they influence total energy sharing and consequently (3- electrons energy spectrum. 
The (3- electrons differential energy spectrum can be written in the following form 0: 

: AF(pfs, Z)E m £ n P n ( P p)e n (e 2 n - rnyf 2 9(e n - m„c 2 ), 

(2) 

e n = W - e p - E n , ep = Ep- m e c 2 . 



dN(e ) 



The sum in (Q) runs over final states of T 3 He + molecule; Ep, ep are the total relativistic- 
and kinetic energies and pp is the momentum of /3-electrons; F(pp, Z) is the so called Fermi 
factor which takes into account Coulomb charge Z of the daughter nucleus 0]. If we will 
count final energies E n from the ground E g state of T 3 He + , the end-point (3 -electron kinetic 
energy Wq is given by the relation: 

Wo = (M t - M a )c 2 + E T2 — Eg — Er — m e c 2 , (3) 

where Et 2 is the ground state of parent T<i molecule; Er = is the molecular center 
of mass recoil energy. The final states population probabilities P n (pp) are given by the 
corresponding matrix elements: 

P n (pp)=J^\<n\e^\T 2 >\ 2 , q = &, (4) 

where R is the internuclear distance ||. The f3- spectrometers currently operating in Mainz 
and Troitsk detect actually integral spectrum for which we have : 

N(ep) = jF(pp, Z)Eppp En Pn(pp) (e 2 n - m^c 4 ) 3/2 9(e n - m v c 2 ). (5) 

The standard procedure consists of fitting experimental spectrum close to the end-point 
energy Wo using formula (f|) varying three parameters A, E Q , w? v and value of background. 

Since first quantitative results have been published in the literature, immense work has 
been done to purify theoretical spectrum in many different aspects ||. They showed that 
about 99% of final states spectrum covers over 100 eV. It mostly consists of the several eV 
narrow rovibrational multiplets wired on bound and resonance adiabatic electronic states as 
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it is shown in Fig. [|. Some 57% falls on to electronic ground state of T 3 He + . The rovibra- 
tional multiplet consists of hundreds bound and resonance states distributed "...unevenly in 
a somewhat erratic manner" with average rotational quantum number J ~ 22 25 || . The 
reason is a huge, on the atomic scale, recoil momentum q ~ 18 releasing. 

The results of extensive joint fits between experiments and available theoretical tables 
show with a certainty that parameter m 2 v should be set at a negative value. The definite 
value explicitly depends on the chosen energy fit interval. Though physically unsatisfactory, 
it permits to draw up an upper estimate m v < 2 3 eV/c 2 . 

The critical values oim u ~ 1 eV/c 2 will be decisive, e.g. for "dark matter" problem. The 
series of experiments aiming to hit this border are on the list of Mainz- Troitsk collaboration. 
The "negative m^" problem will represent then one of the serious bottlenecks. To exclude 
reasons related to the FSS we went to reanalyze basic issues in the theoretical final states 
spectrum and its influence on the (3 -electrons spectrum. 



As Eqs.(0), ([5]) suggest, the sign of non-zero neutrino mass is the non-analytic behavior 
close to the branch energy of given channel. We would like to show, however, that a some- 
what simpler functional form is enough for all practical purposes. Let us take for granted 
that ml is a small parameter: 



II. /3-SPECTRUM: GENERAL FEATURES. 




(6) 



which gives for the spectral sum itself: 




(7) 



We can further expand ^-function 




(8) 
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to show that neutrino mass term enters into the spectrum, essentially, via linear energy 
term (|7|). We plot in Fig. |2] absolute difference between spectral sum in (|5[) and its linear 
approximation (0) using FSS from 0. It shows that absolute difference is uniformly bounded 
and decreases roughly as (m^c 2 ) 4 / {Wo — ep). The formula (|7|) can be further transformed to 
the particularly transparent form: 



e 3 - 3 < E n > t e 2 + 3 < E* > t e - §m 2 c 4 (e- < E n > e )- < E\ > e \ , (9) 
where we introduced FSS cumulative momenta 

P e = ]T#(e n ), <E n>€ =l/P e J2E n 6(e n ) ... 

n n 

It shows that /3-spectrum can be parametrised using few statistical characteristics of the full 
FSS. 



III. /3-SPECTRUM: OPERATOR FORMULATION. 

To find some way to calculate cumulative momenta directly, we use in (|2|) well-known 
substitution E n — > H, where H is Hamiltonian of T 3 He + molecule. Completing the sum 
over final spectrum, we get: 

3/2 



F{p ,Z)EpppJ^< T 2 |e-^> R ) e 2 - nrr 



2 A 



6(W -ep-H)e i ^\T 2 >. j 



(10) 



here e = W — ep — H. It is easy to justify the following commutation relation: 



(11) 



2M M 

using explicit form of kinetic energy operator Tk = — ^A. We can eliminate now recoil 
exponent, making simultaneous substitution in ([T0|): 



W - ep - H ->• W - ep 



<t 



lq, prj 

2M M 



H 



(12) 
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The term ^ can be interpreted as an integral rotational recoil energy shift. It adds up to 
the center-of-mass recoil energy 



Er — > E R + 



n 



2m e c 2 



(q.PR) 

M 



2M 2M t M t 

The additional operator term in (|12|) is estimated at 
expansion parameter. 

To start with ^-function in (flO|) , we use following relation: 



f(A + e) = f(A) + f'(A)e + -f"(A) [e, A] + 0(e 2 , [A, [e, A]}, 



to get 



(13) 



10 . We will use it like an 
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AI 



(14) 



The general functions must be understood as usually in terms of convolution with smooth- 
class functions. In our case it is naturally provided by the final energy resolution. The 
experimentally observed spectra are defined by 



N exp {e p ) = J deR(ep-e)N(e), 



where N(e) is the one from @. Taking for example generic case of Gaussian 

1 



R(e) 



-e 2 /2a 2 



V2 



TIG 



we conclude that 5-like terms in (|14"D will produce local corrections the type shown in Fig.|2| 
Moreover, due to additional / dq integration leading order correction will be quadratic in 
^ ■ Their total spectral power should be sufficiently small, therefor we keep in (|H|) 
leading term only. 

It is easy to follow that expansion of "3/2" -function over parameter with conse- 

quent J dq average will give: 

(e -m v c) - 



M 



(15) 



The particular form of C depends on the operator order in ( fL5|) and is given by 



c = - 1 -(±\ 

3 \MJ 



d_ 
dR 



d 

dR 



d 

dR 



d_ 
dR 



The spectral power of this term is relative small, e.g. for T 3 He + electronic ground state we 
have \C\ < 0.1 eV 3 . 

Putting all together we get a new representation for the /3-spectrum: 

N(e p ) = Z)E m < T 2 |[(e 2 - m^f 2 - (fjf A R e + 6}6(e)\T 2 > . (16) 

The differential spectrum is written down analogously: 



den 



AF(p p , Z)E m < T 2 \[e (e 2 - m 2 c 4 ) 1/2 - | (^f A R ]9(e)\T 2 



> 



(17) 



The new final states expansion formula appears after plug-in J2n \n >< n\ = 1 into ( 161) 
(|17|). For the integral spectrum we have: 

N(e p ) = jF( Pf3 ,Z)E m EnPn[(e 2 n -mlc 4 f 2 - (§f e n < T 2 \A R \n >< n\T 2 > / P n + 



< T 2 \C\n >< n\T 2 > /P n ]9(e 



n ) i * n 



P, 



< To\n > 



(18) 



IV. DISCUSSION 

Let us discuss derived formulas. Compare formulas (^) and ( |TED we get the following 
operator expressions for the cumulative momenta: 

P e =<T 2 \6{e)\T 2 > , <E n>e =<T 2 \H6{e)\T 2 > /P e + ^ 

2 ( / 

< El > e = [< T 2 \(H + £f) 2 8(e)\T 2 > - (&) < T 2 \A R \T 2 >}/P e . 
As it follows from (pT8|) , the ^settc/o-spectrum of final excitations contains vibrational 

2 

modes only. The term represents uniform rotational recoil energy shift. Consider as 
an example T 3 He + ground electronic state. The pseudo-spectrum of final excitations is 
emptied by the overlap with the first few vibrational states 
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I < T 2 \g,v = > | 2 = 52.2%, | < T 2 \g,v = 1 > | 2 = 4.62% 
| < T 2 \g, v = 2 > | 2 = 0.45%, | < T 2 \g, v = 3 > | 2 = 0.08% . . . 

The total sum, as before, is equal to the integral probability 57.4%. Using calculated pop- 
ulation probabilities, we can find electronic ground state first energy moment 

< En >g= 2M + V^E^n = 1-72 + 0.03 = 1.75 eV. 

To get exactly same result one has to sum some hundreds lines in the standard approach. 
Note that vibrational recoil shift is generally ~ 0.1 eV. The sign of the shift depends on the 
particular state. As for the first excited electronic state it is negative. In the relation to the 
highly excited electronic states our approach also give clear recipe on how to incorporate 
correctly nuclear motion. It will be enough to consider that in the final state nuclei move 
under effective Coulomb repulsion Zizz 2)/R. 

We would like to discuss now one very interesting consequence of our formulas. We have 
noted already in the course our derivation that composite center-of-mass, rotational recoil 
shift is 

R M t \ 2m e c 2 ) M t 
That corresponds to the effective change of the end-point energy W as it is seen at W — 
below true end-point energy W : 

wt ff) = W -5W , SW = ^p-. 

M t 

One can guess that Wo fit value will be dominated by the lower energy end, where spectral 
yield is huge. At ~ 200 eV below end-point energy we have 5W ~ 0.04 eV. This value 
is comparable with the statistical uncertainty of W . Fixing Wq & ^\ we sweep into the W 
energy region, where small neutrino mass term is comparable with the total spectral yield. 
To compensate deficit of spectral yield caused by lower value of , we have to make m 2 

negative. The rough estimate give us some few eV 2 of negative shift. The fit experiments 
have to show real size of the effect. 



The integral recoil energy shift is clearly associated with the p/3-dependence of the fi- 
nal states probability distribution P n {pp). The precise nature of that dependence wasn't 
available, apart the fact that P n (pp) change wildly at small variations of q{pp)- 

Last notes concern of the other possible experimental situations. The case of TH de- 
caying molecule passes through our consideration with obvious changes. To note some: the 



rotational recoil shift 2 M t (i+Mt/M ) * s a bout half of the T 2 case; the composite recoil shift 
(|13|) stays the same. The case of initial nonzero J rotational state leads in final channel to 
the selection of J rotational states, which causes uniform spectral shift ~ — 2Mb? — ' Finally; 
there are shouldn't principal difficulties to make proper modifications caused by the solid 
state surrounding. 

As a conclusion, we can suggest one modification which is easy to incorporate into 
standard spectrum. It is slight energy dependence of the end-point energy 

This modification is almost irrelevant for the Wq itself but it can sensibly change final value 
of ml parameter. 
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FIGURES 



FIG. 1. The final states spectrum. 

FIG. 2. The absolut difference between spectral sum and its linear approximation for m v = 1 eV 
reference value. The additional line shows (m u c 2 ) 4 /(Wo — ep) trend line. 
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